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Abstract. Recent mass measurements show a substantial weakening of the binding energy difference 
5 2p {Z, N) = E(Z - 2, N) - 2E(Z, N) + E(Z + 2, N) in the neutron-deficient Pb isotopes. As S 2p is often 
attributed to the size of the proton magic gap, it might be speculated that reduction in 5 2p is related to a 
weakening of the spherical Z = 82 shell. We demonstrate that the observed trend is described quantitatively 
by self-consistent mean-field models in terms of deformed ground states of Hg and Po isotopes. 

PACS. 21.10.Dr - 21.60.Jz - 21.10.Gv - 27.70. +q - 27.80. +w 
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1 Introduction 

The weakening or "quenching" of spherical shell closures 
when going away from the valley of stability to weakly 
bound nuclei is a phenomenon of great current interest. 
It is now well established for neutron-rich N = 20 and 
N = 28 isotones. The study of excited states of nuclei 
around 32 Mg populated by P~ -decay Q strongly indicates 
that neutron-rich N = 20 nuclides are deformed. This is 
consistent with direct mass measurements which, at magic 
numbers 20 and 28, do not yield the familiar drop in 
the two-neutron separation energy (^,|||§]. Energies and 
B(E2) transition probabilities obtained from Coulomb ex- 
citation of radioactive beams of nuclides around 44 S demon- 
strate also that N — 28 isotones below 48 Ca show strong 
collectivity in contradiction with a spherical shell closure 
|5j . Also, there are experimental hints from f}~ decay stud- 
ies of 80 Zn H and the systematics of 2 + and 4 + excitation 
energies in Cd and Pd isotopes that the N = 50 and 
N = 82 shells are weakened when going towards neutron- 
rich nuclei. Shell-quenching in neutron-rich systems has 
far-reaching consequences for astrophysics as it influences 
the r-process path ||. 

The actual sequence of magic numbers in neutron-rich 
nuclei is strongly influenced by an increasing diffuseness 
of the neutron density, the closeness of the particle contin- 
uum, and/or the changes in the spin-orbit splitting P,|lO||. 

a Present address: Service de Physique Nucleaire Theorique 
et Physique Mathematique, CP229, Universite Libre de Brux- 
elles, B-1050 Bruxelles, Belgium. 



For light neutron-rich N = 28 isotones, this leads to an in- 
creased collectivity and even strong stable quadrupole de- 
formation; see, e.g., refs. jll 12,lj|. (For more discussion 
of shell structure of neutron-rich magic and semi-magic 
nuclei, the reader is referred to ref. [ fl4| where a more 
complete literature review is given.) 

The experimental signatures of "magicity" could some- 
times be contradicting. Consider, e.g., the N = 40 sub- 
shell in 68 Ni. While it is clearly visible in the systematics 
of excited states |u||l6|], no sign of a shell effect is found 
when looking at two-neutron separation energies. Theo- 
retically, this apparent inconsistency can be explained in 
terms of dynamic correlations beyond the mean field @. 
The calculated single-particle spectra in 68 Ni indeed show 
a N = 40 shell. However, the sudden increase in collec- 
tivity in the adjacent isotopes obscures the picture when 
looking at isotopic energy differences because one now 
compares nuclei with different intrinsic structure. This ex- 
ample demonstrates that the various signatures of shell 
closures are not always equivalent. 

All examples mentioned so far concern neutron shells. 
The situation seems to be different for protons. For light 
nuclei there is no indication that the proton shell closures 
fade away towards the proton drip line. The analysis of 
recent large-scale mass measurements of proton-rich nuclei 
around Z = 82 at GSI |51| indicates a weakening of the 
proton shell gap in neutron-deficient Pb isotopes. This 
is found when looking at the binding-energy differences 
[191 and Q a values |2fJ. It is the aim of this paper to 
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analyze these data from a theoretical perspective using 
self-consistent mean-field models. 



2 Signatures of shell closures 

The notion of a shell closure comes from a mean-field de- 
scription where one has full insight into the single-nucleon 
energies as eigenvalues of the single-particle Hamilto- 
nian. A shell closure is associated with a large gap in the 
spectrum of e&. It has to be stressed that Ck are not equiva- 
lent to experimental single-particle energies Sk which are 
measured as one-nucleon separation energies. To calcu- 
late those, one has to take into account residual inter- 
actions (such as pairing and coupling to low-lying vibra- 
tional modes), as well as rearrangement effects and the 
core polarization due to the unpaired nuclcon; sec, e.g., 
refs. H|||. 

Bunchiness of levels and the presence of gaps in the 
single-particle spectrum can be quantified in terms of the 
shell-correction energy 



-Bshell = 



- J deeg{e) 



(1) 



where g(e) is a smoothed level density. -E s hcii is large and 
negative for magic nuclei because there the discrete sum 
has a particularly low value while the smoothed expression 
varies only slowly with mass number. It is to be noted that 
-E'shoii is a measure of the deviation of the actual level den- 
sity at the Fermi energy from the smoothed level density. 
For very heavy systems beyond lead, the large level den- 
sity inhibits pronounced shell gaps, and, yet, a dilution of 
levels suffices to produce shell stabilization without magic 
gaps; see ref. Shell corrections discussed in this study 
were calculated using the Green's function approach of 
ref. |2j| . The advantage of this new method is the proper 
treatment of unbound states which appear close to the 
Fermi level in weakly bound (proton-rich or neutron-rich) 
nuclei. 

Neither nor -E s heii are directly measurable quanti- 
ties. A quantity which is accessible from mass systematics 
is the so-called two-proton shell gap 

S 2p (Z, N) = E(Z - 2, N) - 2E(Z, N) + E{Z + 2, N) , 



which is the approximation to the second derivative of 
the nuclear binding energy. Let us assume that the single- 
particle energies do not change within the three isotones 
considered in Eq. (Q), and that the change in the total 
binding energy comes from the variation of occupations 
around the Fermi surface. In such a case, Koopman's the- 
orem states that 52 P should represent twice the gap in the 
corresponding single-particle spectrum. This requires, how- 
ever, that no dramatic rearrangements happen among the 
three nuclei involved in 52 P - For many nuclei, such rear- 
rangements are small and S2 P is a good representation of 
a shell gap. But this interpretation does not hold in situa- 
tions where adding or removing just two nucleons induces 



a substantial change of the mean field. In spite of this 
weakness, it is worth noting that binding-energy differ- 
ences involving even-even nuclei only - like £>2 P and, to 
some extent, Q a - provide a cleaner signal and spectral 
properties than the binding energy expressions connect- 
ing to the neighboring odd-mass systems. The double dif- 
ference 5 P with odd-mass neighbors mixes shell gap, and 
possibly rearrangement effects, with pairing effects, which 
surely imply a dramatic rearrangement of the pairing field 
by blocking. 

Other experimental signatures of magic numbers are 
excitation energies of vibrational states and the associ- 
ated electric transition rates. The gap in the single-particle 
spectrum sets the scale for the lowest excitations. Peaks 
in the systematics of the lowest 2 + states in even-even 
nuclei reflect the stiffness of the potential energy surface 
which is large in closed-shell nuclei. This signature is some- 
what masked by the residual interaction which shifts the 
excitation energy considerably below the lowest particle- 
hole (or two-quasiparticle) energy. Compared to binding- 
energy differences, the data on vibrational states have the 
advantage that they do not mix information from differ- 
ent nuclei. On the other hand, calculations of collective 
excitations require the use of techniques which go beyond 
the mean-field approach; hence the information on magic 
gaps is very difficult to extract. 



3 Theoretical framework 

We investigate the stability of the Z = 82 shell in the 
framework of the self-consistent mean-field theory, cither 
non-relativistic, using Skyrme interactions (SHF), or the 
relativistic mean-field (RMF) approach. As earlier stud- 
ies have shown that paramctcrizations of the models of- 
ten differ when extrapolated to weakly bound systems, 
we choose a sample of representative parameterizations 
which all give a very satisfactory description of stable nu- 
clei but differ in details. Namely, we used the Skyrme 
interactions SkP @, SLy6 @, SkI3, and SkI4 jp. In 
the relativistic calculations, we employed NL3 [p9| and 
NL-Z2 |50) parameterizations. The force SkP has effective 
mass to* /to = 1 and was originally designed to describe 
the particle-hole and particle-particle channel of the ef- 
fective interaction simultaneously. The forces SLy6, SkI3, 
and SkI4 stem from recent fits which already include data 
on exotic nuclei. Both SkP and SLy6 use the standard 
spin-orbit interaction. The forces SkI3/4 employ a spin- 
orbit force with modified isovector dependence. SkI3 con- 
tains a fixed isovector part analogous to the RMF, whereas 
SkI4 is adjusted allowing free variation of the isovector 
spin-orbit term. The RMF force NL-Z2 is fitted in the 
same way as SkI3 and SkI4 to a similar set of observables. 
Pairing is treated within the BCS approximation using a 
zero-range delta force with the strength adjusted for each 
mean- field parameterization as described as in ref. |3l)|. 
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Fig. 1. Two-proton gap parameter &2 P (|2|) for the chain of 
Pb isotopes obtained in several spherical mean-field models. 
Experimental values are marked with filled diamonds. 
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Fig. 3. Spherical proton shell correction for the chain of Pb 
isotopes extracted from self-consistent calculations. 



4 Results 

4.1 Spherical calculations 

We start with a presentation of purely spherical calcula- 
tions. By constraining the shape to be spherical, one can 
dramatically reduce polarization effects. Consequently, the 
assumptions behind the Koopman's theorem can be met, 
and the clear correspondence between (5 2j , and single-particle 
energies emerges. Figure^ shows the two-proton shell gap 
parameter along the chain of Pb isotopes. The experi- 
mental data are compared with results of spherical mean- 
held calculations. While the experimental values decrease 
monotonously when going from 208 Pb towards the proton- 
drip line, theoretical results show a very different trend. 
Namely, S2 P slightly decreases with Z. The corresponding 
spherical single-proton energies in SkI3 are displayed in 
fig. 0. It is seen that the Z = 82 shell gap stays large for all 



2f7/2-..; nj ' 

_ lh9/2~,, \. n/""\'*x* 


N N$^ SkI3 1 




.,'••...'*••-., lill/2 + 




-•N. api/l 


3sl/2 

2d3/2 + ?\ 
■lhll/2 

2d5/2 + ^ 
lg7/2 + ^ 


*\/\ \ "'- 3 P :i 2 
\\. "'"-2f5/2" 

' • '" •V s lil3/2 ' 
2f7/2 

\ ® , " 
\. --lh9/2 


.lg9/2 + -^ 


^\ 3sl/2 + 
\^ ' ^Sn 2d3/2 

^^2d5/2 + 


2pl/2"-/ S \^ 



I I I I I I I 

90 100 110 120 130 140 
Neutron Number N 



Fig. 2. Self-consistent single-proton energies ek for the chain 
of Pb isotopes calculated with the spherical SkI3 model. Other 
effective interactions give similar results. See ref. |}o| for a 
comparison of calculated and experimental spectra for 08 Pb. 



the Pb isotopes considered, and it grows slightly when ap- 
proaching the proton drip line. (Other parameterizations 
employed in our study give similar results.) This result is 
consistent with fig. ffl: S2 P reproduces twice the shell gap 
seen in the spectrum of This demonstrates that the 
idea beyond the two-nucleon shell gap as a signature of 
shell closures is correct, provided the mean field does not 
change much between the nuclei involved in its evaluation. 
The restriction to spherical calculations enforces this fea- 
ture. 

The stability of the Z — 82 shell is confirmed when 
inspecting the proton shell correction energy -Esheii.p ex- 
tracted from our self-consistent calculations; see fig. |j| 
There are some differences among the various forces as 
far as the actual values are concerned. For instance, RMF 
models generally yield smaller shell correction than SHF 
models. The common feature seen in fig. || is that all the 
models predict the same trend, namely that the magnitude 
of the proton shell correction is slightly increasing when 
going proton-rich. One can thus say that all the indicators 
of shell effects considered in our study, i.e., single-particle 
energies, shell corrections, and 52 P , give a consistent pic- 
ture of a stable Z = 82 shell for all neutron-deficient Pb 
isotopes. Contrary to what may be suggested by the ex- 
perimental plot of (52 P , the proton magic gap does not de- 
crease, but slightly increases when approaching the proton 
drip line. 

4.2 Deformed calculations 

It is well known that shape coexistence is a common fea- 
ture in neutron-deficient Hg, Pb, and Po nuclei, and that 
most of the Hg and Po nuclei have deformed ground states 
[ ^2p3| . Figure |I] shows the calculated deformation energies 
AE (defined as a difference between ground-state bind- 
ing energies obtained in deformed and spherical calcula- 
tions) and quadrupole deformations for even-even Hg, 
Pb, and Po isotopic chains. Most nuclei have deformed 
ground states with an energy gain of several MeV. The 
potential energy surfaces calculated in this region exhibit 
several competing minima, often leading to shape coex- 
istence and configuration mixing, see e.g., refs. ^J^l . 
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Fig. 4. Calculated deformation energy AE (top) and associated quadrupole deformation = 4ty (r 2 Y20) / {3Ar 2 ) with 
ro = 1.2 A 1 ' 3 fm (bottom) for the ground states of even-even Hg, Pb, and Po isotopes. Experimental quadrupole deformations 
(filled diamonds) are taken from ref. f49]. 



In fact, the question if whether there were strongly de- 
formed nuclei near the magic proton shell Z = 82 was al- 
ready raised as early as 1972 |37],|38|]. Experimentally, the 
even-even 176 -i90jjg i so t p es have oblate ground states 
but also prolate excited states |3(|. The data on excita- 
tion spectra and charge radii for the Pb isotopes, on the 
other hand, are consistent with the assumption that the 
ground states down to N — 90 are spherical [^3|j39| , ^0|| , al- 
though one has to be aware that the potential landscape 
for Pb also becomes softer for proton-rich isotopes. As 
a result, coexisting oblate and prolate minima show up 
|32|,[35|, the most spectacular being the multitude of co- 
existing structures in 186 Pb A systematic analysis of 
available excitation data in Po isotopes demonstrates that 
collectivity increases rapidly in neutron-deficient isotopes 

HI)- 

When deformation effects are taken into account, the 
behavior of b~2p changes dramatically. As shown in fig. 0, 
all our models predict gradually decreasing 52 P towards 
the proton drip line in accordance with experiment. The 
crucial point is that one deals with systems where the 
mean field is extremely sensitive to small changes in nu- 
clcon number. These rearrangement effects strongly affect 
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Fig. 5. Two-proton gap parameter 52 P (Q) for the chain of 
Pb isotopes obtained in several deformed mean-field models. 
Experimental values are marked with filled diamonds. 



the behavior of 82^ and mask the presence of the magic 
Z = 82 proton shell; it is nuclear deformation that ex- 
plains the apparent discrepancy seen in fig. [l| 

Complementary information on shell effects can be drawn 
from a-decay studies. An analysis of reduced a widths and 
Q a values of nuclei around 182 Pb suggests a weakening of 
the Z = 82 shell effect [Q. However, since Q a values are 
obtained from a finite-difference formula similar to S2 P , 
experimental systematics Q a reflect deformation effects 
and the structural change of ground states between the 
isotoncs of Hg, Pb, and Po. Moreover, a-decay hindrance 
factor systematics can be understood assuming the stabil- 
ity of the spherical Z = 82 shell at the very proton-rich 
side |43|,H and can be described in the language of shape 
mixing [03Eq]. 

The rich structure of potential landscapes in the region 
of proton-rich Pb nuclei provides a challenging benchmark 
for theoretical modeling. The quantitative description of 
all details is beyond the abilities of current self-consistent 
mean-field models and, because of the presence of low- 
lying intruder states, also requires the use of sophisti- 
cated techniques which incorporate configuration mixing 
effects. Small shifts in single-particle levels or changes in 
the parameterization of pairing interaction do influence 
relative positions of coexisting minima. Therefore, many 
forces miss the exact location of the onset of deforma- 
tion for the Hg isotopes, see fig. ^ and [Q. Some forces 
(RMF and SkI3) even predict deformed ground states in 
some Pb isotopes - in contradiction with experiment. We 
have seen in fig. ^ that the inclusion of deformation im- 
proves the agreement of theory with experiment for #2 P - 
Still, some discrepancies remain. Note that all SHF forces 
already overestimate #2p for the doubly magic 208 Pb, al- 
though this nucleus is spherical, as well as its neighboring 
nucleus. The simplest excuse is that the single-particle 
states in 208 Pb are not perfectly described by the mod- 
els (see, e.g., p0|). However, another important source of 
discrepancy are correlations beyond the mean field. Their 
influence on the two-neutron separation energies around 
208 Pb has been studied using the generator coordinate 
method (GCM) in ref. JitJ. The important outcome is 
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that the quadrupole correlations decrease the two-neutron 
shell gap at the shell closure. A similar effect can be ex- 
pected for the two-proton shell gap parameter. The effect 
of quadrupole correlations on 52 P is probably largest in 
the transitional region of the onset of static deformation 
around N = 120 where also the discrepancy between ex- 
periment and our mean-field calculations is largest. Con- 
figuration mixing effects in proton-rich Pb isotopes were 
studied in ref. |4q| , but were aimed at excitation energies, 
and not mass systematics (see also recent GCM work ]50|). 



5 Summary 

Finite-difference binding-energy indicators such as &2 q or 
Q a lose their validity as signatures of shell closures as soon 
as they are used in a region where the structure of nuclear 
ground states is rapidly changing. Such changes of nuclear 
ground state configurations with Z and N explain quali- 
tatively the systematics of experimental S2 P and Q a in the 
neutron-deficient Pb region. In particular, our calculations 
do not support previous speculations about quenching of 
the Z — 82 shell in the neutron-deficient Pb isotopes. 

It is important, however, to remember that the com- 
plex structure of nuclei from shape-coexisting regions lim- 
its their description in terms of mean-field models. The 
presence of low- lying intruder states having different shapes 
requires the use of theories which can account for config- 
uration mixing effects. Therefore, the whole notion of a 
"shell gap", which originates from the mean-field theory, 
can be questionable in such situations. In general, one has 
to be aware that the usual signature of spherical shell clo- 
sures might not be robust when going to exotic nuclei far 
from stability in which correlations play an important role. 
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